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ABSTRACT 

A growing number of observations reveal a subset of Type la supernovae undergoing circumstellar 
interaction (SNe la-CSM). We present unpublished archival Spitzer Space Telescope data on SNe 
la-CSM 2002ic and 2005gj obtained > 1300 and 500 days post-discovery, respectively. Both SNe 
show evidence for late-time mid- infrared (mid-IR) emission from warm dust. The dust parameters 
are most consistent with a pre-existing dust shell that lies beyond the forward-shock radius, most 
likely radiatively heated by optical and X-ray emission continuously generated by late-time CSM 
interaction. In the case of SN 2005gj, the mid-IR luminosity more than doubles after 1 year 
post-discovery. While we are not aware of any late-time optical-wavelength observations at these 
epochs, we attribute this rebrightening to renewed shock interaction with a dense circumstellar shell. 

Subject headings: circumstellar matter — supernovae: general — supernovae: individual (SN 2002ic, 
SN 2005gj) — dust, extinction — infrared: stars 



1. INTRODUCTION 

The ability to standardize Type la supernova (SN la) 
light curv es yields preci se cosmological distance indica- 
tors fe.g.. IPhillipslll993[ ). Despite these empirical rela- 
tionships, questions remain about the underlying physics 
and progenitor systems. The SN itself is generally ac- 
cepted to be the thermonuclear explosion of a C/0 white 
dwarf (WD), but the nature of the companion star re- 
mains ambiguous. Evidence now exists for both single- 
dege nerate and double-degenerate (SD, DD) channels 
Ce.g. iPatat et al.ll2007l: iBloiidin et a"l]l2009l: ISmon et alJ 



2009i: iNugent et all 1201 iMGaneshalingainet all 120111: 



Brown et all l2012t iFolev et all 12012 
2012HSilverman et al.ll20ll 



ggam et ai.l l/uiit 
alibi: iBloom et al.l 



Recent studies reveal a subsample of SNe la that ex- 
hibit signs of significant interaction between the for- 
ward shock and a dense circumstellar me dium (here- 
after SNe la-CSM; ISilverman etHI I2013U and refer- 
ences within). The spectra show many similarities to 
those of SNe Iln, defined by their strong and relatively 
narrow hydrogen emissio n lines generat ed fr om a slow, 
dense, pre-existing CSM (|Schlegell 119901 : see iFilippenkol 
119971 for a review). The SNe la-CSM also have slightly 
larger peak luminosities than typical SNe la, in the range 
—21.3 < Mj^ < —19 mag. Detailed observations of 
the nearby SN la-CSM PTFllkx confirm a probable SD 
channel for at least this par ticular event (|Dildav et al.l 
I20TI ISihiFman et al.ll20T3al) . 

Two of the m ost well-studied SNe la-CSM are 
2002ic and 2005g.i (iHamuv et all 120031: iDCTTg et al.ll2004l: 
Kotak et al]|200l IW ang et al."2004': 'Wood-Vasev et al.l 
2004; Ald cring et al . 2006; Prieto et al. 2007). The data 
suggest significant mass loss (of order > 10^** M© yr^^) 
from the companion star and significant amounts of 
warm dust emitting in the near-infrared (near-IR) . These 
characteristics are also similar to those of SNe Iln. Ow- 
ing to their dense CSM, SNe Iln exhibit late-time (> 100 
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day) IR emission from warm dust that is continuously 
heated by visible and X-r ay radiation generated by ongo- 
ing C SM interaction (e-g.- ICerardv et al.ll200llFox et al] 
I2OIII ). The CSM geometry derived from these dust shells 
reveals important clues about the progenitor mass-loss 
history. 

In this Letter we present unpublished archival 
Spitzer Space Telescope data on SNe 2002ic and 2005gj 
obtained > 1300 and 500 days post-discovery, respec- 
tively. Section |2| lists the details of the observations; 
Spitzer photometry constrains the dust mass and tem- 
perature, and thus the luminosity. We explore the origin 
and heating mechanism of the dust in fj3l fj4| presents our 
conclusions. 



2. OBSERVATIONS 



2.1. 



Warm Spitzer /IR AC Photometry 

A search in the Spitzer Heritage Archive (SHAJ^ 
revealed unpublished observations of SNe 2002ic and 
2005gj , summarized in Table [T] The SHA provide ac- 
cess to the Post Basic Calibrated Data (pbcd), which 
are already fully coadded and calibrated. Figure [T] shows 
false-color images of the combined 3.6, 4.5, and 5.8 /im 
channel images at a single epoch. The background flux in 
most of the SN host galaxies is bright and exhibits rapid 
spatial variations. Although template subtraction is a 
commonly used technique to minimize photometric con- 
fusion from the underlying galaxy, no pre-SN Spitzer ob- 
servations exist. Instead, photometry was performed us- 
ing the DAOPHOT aperture-photometry package (APPHOT) 
in IRAfB SN 2005g.i falls < 1" from the galactic nucleus 
(|Aldering et all 120061 : iPrieto et all I2007D . To minimize 
contributions from the underlying galaxy in both cases. 



SHA 



be 



accessed 



from 
http://sha.ipac.caltech.edu/applications/Spitzer/SHA/ . 

* IRAF: the Image Reduction and Analysis Facility is distributed 
by the National Optical Astronomy Observatory, wrhich is oper- 
ated by the Association of Universities for Research in Astronomy 
(AURA) under cooperative agreement with the National Science 
Foundation (NSF). 



Table 1 

Spitzer Observations^ 



SN 


JD 


Epoch 


PID 


AOR 


Distance 


tint 


3.6 ^m^ 


4.5 /im'' 


5.8 /xm^ 


8.0 ^m^ 




-2,450,000 


(days) 






(Mpc) 


(s) 




(lOi'^ erg S-- 


' cm-2 A-i) 




2002ic 


3386 


795 


3248 


10550272 


280 


4000 


0.150(0.075) 


0.136(0.057) 


0.092(0.038) 


0.052(0.023) 


2002ic 


3770 


1179 


20256 


14455040 


280 


3600 


0.060(0.047) 


0.063(0.039) 


0.052(0.029) 


0.040(0.020) 


2002ic 


3961 


1370 


30292 


17965824 


280 


2400 


0.042(0.039) 


0.045(0.033) 


0.037(0.024) 


0.031(0.017) 


2002ic 


4356 


1765 


40619 


23107840 


280 


2400 


<0.03 


<0.03 


<0.03 


<0.03 


2005gj 


3778 


139 


264 


16868096 


268 


2400 


0.130(0.069) 


0.092(0.047) 


0.036(0.024) 


0.017(0.013) 


2005gj 


4004 


365 


30733 


19308800 


268 


2400 


0.183(0.082) 


0.110(0.052) 


0.049(0.028) 


0.024(0.015) 


2005gj 


4149 


510 


30733 


19309056 


268 


2400 


0.276(0.101) 


0.199(0.069) 


0.113(0.043) 


0.049(0.022) 



^ Upper limits for nondetections were derived by the point-source sensitivity in Table 2.10 of the IRAC Instrument Handbook, 

version 2. 

2 All distances are derived from the host-galaxy rcdshift assuming Hq = 72 km s^-"^ Mpc^^. 

^ la uncertainties are given in parentheses. 



a 2-pixcl radius was chosen and aperture corrections were 
applied, although removal of nuclear contributions is ad- 
mittedly difficult. 

2.2. Dust Temperature and Mass 

Assuming only thermal emission, the mid-IR photome- 
try provides a strong constraint on the dust temperature 
and mass, and thus the IR luminosity. For optically thin 
dust with mass M^ and particle radius a, at a distance 
d from the observer, thermally emitting at a single equi- 
librium teinperat ure Td, the flux can be written as (e.g., 
IHildebrandlll983l ) 



F„ 



MdB,iTd)KAa) 
d^ 



(1) 



where B^(Td) is the Planck blaekbody function and K^(a) 
is the dust absorption coefficient. 

For simple dust populations of a single size composed 
entirely of either silicates or graphite, the IDL MPFIT 
function ()Markwardtll2009l) finds the best fit (see Figure 
[2]) of Equation [1] by varying Md and Td to minimize the 
value of Y^ . The absorption coefficients, k, are given in 
Figure 4 of iFox et all (I2010D . We do not include the 8.0 
/im data in our fits because of an apparent excess (see 
below). With only three data points at each epoch, we 
limit our fits to a single component (see Figure [5]). Table 
[21 lists the best-fit parameters for graphite gra ins of size 
a = 0.1 /im (standard grain parameters; see iFox et alj 

MM- 




2.3. Origin of the 8 ^m Excess 

There are several potential origins for the excess flux 
observed at 8.0 //m in Figure [H (1) The 8.0 /im data 
may be the actual SN dust component, while the shorter 
wavelengths are dominated by an underlying galactic 
component. No single model, however, can reproduce 
the shorter-wavelength behavior in both SNe 2002ic and 
2005gj. We therefore argue that the shorter- wavelength 
data are dominated by SN dust. (2) Spitzer's Channel 
4 bandpass extends out past 9.0 //m, which includes the 
blue wing of a strong 10 /jm feature a ssociated with sil- 
icate dust grains (see IFox et all I2010D . We tried fitting 
silicate dust models to the Spitzer data in Figure [2J but 
find that the 8.0 /xm data are not consistent with the 
expected intensity. We did not account for the Channel 
4 filter transmission curve, but we do not expect this to 
have a large impact since the transmission remains fairly 
constant out past 9. We therefore rule out silicate dust, 
but note that it may contribute if a small fraction of the 
dust consists of silicates. (3) The 8.0 //m data may be 
caused by a second, colder SN dust component peaking at 
longer wavelengths. (4) The excess flux may be a result 
of a second, colder galactic dust component, but difficult 
to remove given the SN positions relative to their respec- 
tive galactic nucleus and the poorer resolution at longer 
wavelengths. Because we lack longer-wavelength data, 
we cannot discriminate between these final two possibil- 
ities. Regardless of the origin, we exclude the 8.0 //m 
data from our fits. 




Figure 1. False-color 3.6, 4.5, 5.8 /im images of SNe 2002ic and 2005gj at late-times. 
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Figure 2. Photometry of SNo 2002ic and 2005gj in Spitzer /IRAC Channels 1 (3.6 fim), 2 (4.5 fim), 3 (5.8 fim), and 4 (8.0 ^m). 
Overplottcd arc the resulting best fits of EguationTll 



3. ANALYSIS AND DISCUSSION 

Figure [3] plots the corresponding IR luminosity evolu- 
tion for each object. SN 2002ic remains bright (L^ > 
10* Lq) for more than 2 yr post-discovery, but continues 
to fade throughout the observations. By contrast, SN 
2005gj brightens after 1 year post-discovery. While we 
only have mid-IR photometry available at these epochs, 
we explore the constraints these data provide on the cir- 
cumstellar environment. 



3.1. Possible Origins and Heating Mechanisms 

The source of the mid-IR emission is warm dust, but 
the origin and heating mechanism of the dust are less 
clear. The dust may be either newly formed or pre- 
existing, andeither shock or radiatively heated; see 
iFox et ail (|201CI| ) for a full discussion. To discriminate 
between possible scenarios, we assume a spherically sym- 
metric, optically thin dust shell and calculate the black- 
body radius, rbb = [Ld/ i4:'KaT^)]^''^ , which sets a min- 
imum shell size. A luminosity Ld ~ 10* L© and dust 
temperature T^ ~ 500 K yield a blackbody radius of 
rbb ~ 10^^ cm. This radius is larger than that of a 
forward shock moving at Vg ~ 10,000 km s^^ on day 
800. This point seemingly rules out, particularly for SN 
2005gj, the possibility of new dust formation in either 
the more slowly moving ejecta or the cold, dense shell 
(CDS) that can form behind the forward shock. 

The alternative to newly formed dust is a pre-existing 
shell. Again, the fact that the blackbody radius is be- 
yond the forward-shock radius rules out the likelihood of 
shock heating. Furthermore, an IR light-echo scenario 
(e.g., [Dwek lOS^), in which the dust shell is heated by 
the peak SN luminosity, is unfeasible. The implied shell 
radii (rocho = ipiatoau/2c) would require peak luminosi- 
ties of nearly 10^^ L© to heat the dust to the observed 
temperatures. 

In fact, the observed dust-shell parameters (i.e., ra- 
dius, tempera ture, mass) are compa rable to those seen in 
SNc Iln (e.g.. iFox et"aLll201lL [20131) . In these cases, the 
pre-existing dust shell is continuously heated by optical 
and X-ray emission generated at inner radii by ongoing 
interaction between the forward shock and dense CSM. 
Assuming an optically thin dust shell, the observed dust 
temperature (Td) and shell radius (rj) require a com- 



bined optical, ultraviolet, and/or X-ray flux given by 



L. 



opt/UV/X 



—par^aTsj^ 



J B^{Tsjq)Qahsii^)dv 



(2) 



for a dust bulk (volume) density p and an effective SN 
blackbody temperature Tsn, where B^ is the Planck 
blackbody function, Qabs is the dust absorption efH- 
cien cy, and k(v) is th e dust absorption coefficient. Figure 
4 in iFox et al.l (|2013[ ) shows that the blackbody radii of 
SNe 2002ic and 2005gj [rbb ~ (0.5-1.0) xlO^^ cm] at 
temperatures Td ~ 500-750 K require optical and/or X- 



ray luminosities in the range 10 < L, 



opt/UV/X 



< 10^ Lr. 



While we don't have optical or X-ray observations at 
these epochs, the most recent optical observations of SNe 
2002ic and 2005gj are consistent, with measured lumi- 
nosities of 10^'^ and lO^"*^ L © on days 25 and 149, re- 
spectively (iDeng et al.ll2004l: iPrieto et al.ll2007l) . 

3.2. Evidence for Shells 

Mid-IR wavelengths probe the characteristics of the 
CSM at the dust-shell radius. Assuming a dust-to-gas 
mass ratio expected in the H-rich envelope of a massive 
star, Zd = Md/Mg w 0.01, the dust-shell mass can be 
tied to the progenitor's total mass-loss rate. 



Mo 



Md 
Zd^r 
3/M, 
4VMo 



^V "^ )( 

I©/ V120kms-i/ V 



5 X lO^*' cm 



(^)Mo yrtri) 



for a progenitor wind speed Ww. The relatively narrow 
lines observed in SNe la-CSM originate in the slow pre- 
shocked CSM and can be used to approximate the pro- 
genitor wind speed. The precursor wind velocities for 
SNe 2002ic and 2005g.i are ^;w = 100 and 60 km s- \ re- 
spectively (jKotak et al.ll2004l: lAldering et al.ll2006D . As- 
suming a thin shell, Ar/r = 1/10, wind speed Uw — 
60 km s~^, and radius rbb = 10^^ cm, the approxi- 
mate mass-loss rate to produce the observed dust shell 
is Afouter ~ 10~^ M© yr~^. A smaller dust-shell radius 
would require an even larger mass-loss rate. 

The optical and/or X-ray emission generated by CSM 
interaction traces the mass loss at the inner radii. As- 
suming a density oc r~^ wind profile, the rate can 



Table 2 

IR Fitting Parameters (a = 0.1 /.tin) 



Days (SN 2002ic) 
800 900 10D0 1100 1200 1300 1400 



Epoch (days) Mj (Mq) T^ (K) Lj (Lq) 



SN 



2002ic 


795 


0.019 


554 


1.54x10* 


2002ic 


1179 


0.019 


489 


7.58x10^ 


2002ic 


1370 


0.013 


488 


5.38x10^ 


2005gj 


139 


0.002 


783 


1.17x10* 


2005gj 


365 


0.002 


805 


1.61x10* 


2005gj 


510 


0.010 


655 


2.35x10* 



be written as a function of the optical/X-ray luminos- 
ity, progenitor win d spee d, and shock v elocity (e.g., 
iChugai fc Danziged 1994; S mith et all 12001 ): 



Mi, 



2v„ 



evi 



-L. 



opt/UV/X: 



L. 



opt/UV/X 



2 1 X 10 "* 

'.3 X 1041 ergs^Wo.5/ 

where e < 1 is the efficiency of converting shock kinetic 
energy into visual light. While the conversion efficiency 
varies greatly depending on shock speed and wind den- 
sity, we assume e ~ 0.5, acknowledging that this value 
may be too high. Again, we do not have late-time op- 
tical or X-ray observations, but we do have theoretical 
estimates from Equation [2] in §3.11 An optical luminos- 
ity iopt/uv/x ~ 10^'^ L©, wind speed v^ = 120 km s~^, 
shock velocity v^ = 10,000 km s"^ (jDeng et al.ll2004[ ). 
and conversion efficiency e = 0.1 correspond to a mass- 
loss rate Minnor ~ 10~^ M0 yr~i. 

While the variables used above are only order-of- 
magnitude approximations, the derived rate for Minnor 
reveals two things about the circumstellar medium. (1) 

The difference between the mass-loss rates (A/jnncr < 
Moutor) suggests that the dust shells were formed dur- 
ing a period of increased, nonsteady mass-loss. (2) Com- 
pared to mass-loss rates derived from optical data at ear- 
lier epochs, our estimate of Mjnner ~ 10""^ Mq yr~i is 
consistent with S N 2002ic on day 250 (for e « 0.1; 
iKotak et al.ll2004 ). but at least an order of magnitude 
larger than that measured for SN 2005gj on day 74 
(jPrieto et al.ll2007l ). The increased circumstellar density 
derived for SN 2005gj at this late time suggests the pres- 
ence of another shell of material. 

The decline of SN 2002ic occurs at >800 days, corre- 
sponding to the time at which the forward shock (vg = 
10,000 km s^i) would reach the blackbody radius (rbb ~ 
10^^ cm). The declining light curve may be attributed 
to the forward shock overtaking and destroying the dust 
and/or a decreasing amount of CSM interaction accom- 
panied by a declining radiative heating source. 

Alternatively, the rebrightening of the pre-existing 
dust in SN 2005gj is likely due to radiative heating by 
renewed shock interaction with this dense circumstellar 
shell. From Equation |4l an order-of- magnitude increase 
in mass loss results in an increase in the optical lumi- 
nosity by a factor of 3-4, assuming the shock velocity 
decreases to 0.8 of its former value. For a constant black- 
body radius of rbb « 5 x 10^^ cm, Equation [2] and Figure 
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Figure 3. The late-time mid-IR luminosity evolution of SNe 
2002ic and 2005gj. The luminosity is derived by summing over 
the modified blackbody fits given in Equation [T] and shown in Fig- 
ure [2] SN 2005gj exhibits an unanticipated rise more than a year 
post-discovery, suggesting renewed shock interaction between the 
forward shock and a pre-existing circumstellar shell. 



10 in lFox et al.l ()2010D show that an increase in the opti- 
cal luminosity from 10* to lO*'* L© results in a dust tem- 
perature increase from '^ 600 to 750 K. Since L^ oc T^, 
this change in temperature results in a luminosity in- 
crease of a factor of 2.4. 

While we observe a luminosity increase of this magni- 
tude (see Tabled , the inferred dust temperature actually 
decreases by ~ 150 K. The implication is that the black- 
body radius must also increase from 5 x 10^^ to 10^^ cm, 
which would explain a lower dust temperature along with 
a higher dust mass and luminosity. The problem with 
invoking this scenario is that the only way to increase 
the blackbody radius (assuming a spherically symmetric 
shell of dust) would be for the increased optical luminos- 
ity to vaporize all dust out to the new blackbody radius. 
For a dust vaporization tem perature o f rcvan ~ 2000 K, 
Equation [2] and Figure 8 in IFox et all ()2010J ) show that 
to vaporize dust out to 10^^ cm requires a luminosity of 
> 10^° L0, which is not likely. The more probable expla- 
nation is that with only three data points, the slopes of 
the curves in Figure [2] and, therefore the temperatures, 
are biased by contamination from the underlying galaxy 
nucleus (< 1" away and a 0.6" pixel scale). The total 
dust luminosity, which is less sensitive to the slope, is 
most consistent with radiative heating by renewed shock 
interaction. 

4. SUMMARY 

This letter presents unpublished archival Spitzer data 
on SNe la-CSM 2002ic and 2005gj obtained > 1300 and 
500 days post-discovery. The mid-IR data show evidence 
of emission from warm dust. While we don't have simul- 
taneous observations at shorter wavelengths, the warm 
dust parameters are most constant with a pre-existing 
dust shell heated by a combination of optical, UV, and 
X-ray emission continuously generated by ongoing CSM 
interaction. The degree of CSM interaction dictates 
the dust temperature and, thereby, the luminosity. In 
the case of SN 2005gj, the mid-IR luminosity nearly 
doubles more than 1 year post-discovery, suggesting an 
increasing amount of CSM interaction. We attribute 
this renewed shock interaction to a dense circumstellar 
shell produced during a period of increased mass loss by 
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the progenitor companion. While progenitor mass-loss 
provides additional support for a SD channel, future 
multi-wavelength observations of SNe la-CSM will be 
necessary to better trace the complete mass-loss history 
and constrain the nature of the companion star. 



This work is based on archival data obtained with 
the Spitzer Space Telescope, which is operated by the 
Jet Propulsion Laboratory, California Institute of Tech- 
nology, under a contract with NASA. Support for this 
work was provided by NASA through an award issued 
by JPL/Cahech (P90031). A.V.F. and his group ac- 
knowledge generous financial assistance from the Richard 
and Rhoda Goldman Fund, the Christopher R. Redlich 
Fund, the TABASCO Foundation, and NSF grant AST- 
1211916. 

REFERENCES 



Aldcring, G., ct al. 2006, ApJ, 650, 510 

Blondin, S., Prieto, J. L., Patat, P., Challis, P., Hicken, M., 

Kirshncr, R. P., Matheson, T., &: Modjaz, M. 2009, ApJ, 693, 

207 
Bloom, J. S., et al. 2012, ApJL, 744, L17 
Brown, P. J., Dawson, K. S., Harris, D. W., Olmstead, M., Milne, 

P., & Roming, P. W. A. 2012, ApJ, 749, 18 
Chugai, N. N., & Danziger, I. J. 1994, MNRAS, 268, 173 
Deng, J., et al. 2004, ApJ, 605, L37 
Dilday, B., et al. 2012, Science, 337, 942 



Dwck, E. 1983, ApJ, 274, 175 

Filippenko, A. V. 1997, ARA&A, 35, 309 

Foley, R. J., et al. 2012a, ApJ, 752, 101 

— . 2012b, ApJ, 744, 38 

Fox, O. D., Chevalier, R. A., Dwek, E., Skrutskie, M. F., 

Sugerman, B. E. K., & Leisenring, J. M. 2010, ApJ, 725, 1768 
Fox, O. D., Filippenko, A. V., Skrutskie, M. F., Silverman, J. M., 

Ganoshalingam, M., Conko, S. B., & Clubb, K. I. 2013, 

^rXiv: 1304.0248 
Fox, O. D., et al. 2011, ApJ, 741, 7 
Ganeshalingam, M., Li, W., & Filippenko, A. V. 2011, MNRAS, 

416, 2607 
Gerardy, C. L., ct al. 2002, ApJ, 575, 1007 
Hamuy, M., et al. 2003, Nature, 424, 651 
Hildebrand, R. H. 1983, QJRAS, 24, 267 
Kotak, R., Meikle, W. P. S., Adamson, A., & Leggett. S. K. 2004, 

MNRAS, 354, L13 
Markwardt, C. B. 2009, in Astronomical Data Analysis Software 

and Systems XVIII, ed. D. A. Bohlendcr, D. Durand, & 

P. Dowlcr (San Francisco: ASP), 251 
Nugent, P. E., et al. 2011, Nature, 480, 344 
Patat, F., et al. 2007, Science, 317, 924 
Phillips, M. M. 1993, ApJ, 413, L105 
Prieto, J. L., et al. 2007, arXiv:0706.4d88' 
Schlegel, E. M. 1990, MNRAS, 244, 269 
Silverman, J. M., ct al. 2012, MNRAS, 425, 1789 
— . 2013a, arXiv:1303.7234 
— . 2013b, arXiv:1304:0763 
Simon, J. D., et al. 2009, ApJ, 702, 1157 
Smith, N., et al. 2009, ApJ, 695, 1334 
Wang, L., Baade, D., Hoflich, P., Wheeler, J. C, Kawabata, K., 

& Nomoto, K. 2004, ApJ, 604, L53 
Wood-Vasey, W. M., Wang, L., & Aldcring, G. 2004, ApJ, 616, 

339 



